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Rotor thrust coefficient 
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Gravity induced body forces per unit mass (m/sec 2 ) 

Impulsive force function; x is the argument. Quantifies impulsive force due to vortex- 
driven flow crossing the slipstream shear layer. 

Number of rotor blades 

Normal vector to control surface, used in vorticity conservation surface integrals 
(positive radially outward) 

Static pressure (N/m 2 ) 

Static pressure component due to quasi-steady vortex-dominated flow (N/m 2 ) 

Static pressure component due to flow with unsteady, non-vortex-driven velocity 
components (dominated instead by large radial and not tangential velocities 
(N/m 2 ) 

Rotor radius (m) 

Radial coordinate, referenced to the vortex core center (m) 

Vortex core radius initial value, at t = 0 (m) 

Outer radial limit of derived stream function solution, where mass transfer across 
slipstream shear layer ceases (m) 


t Time (sec) 

u Velocity component normal to the vortex core control volume surface element 

(positive radially outward) 

U Rotor axial-flight freestream velocity (m/sec) 

Vj Rotor tip speed (m/sec) 

v Mean inflow velocity at rotor disk plane (m/sec) 
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Vi Rotor slipstream velocity in the far-wake, from momentum theory, Vj = 2v (m/sec) 

v r Radial velocity component, with respect to local cylindrical coordinate system with 

origin at vortex core center (m/sec) 

v z Axial velocity component, with respect to local cylindrical coordinate system with 

origin at vortex core center, flow along filament length (m/sec) 

ve Tangential velocity component, with respect to local cylindrical coordinate system 

with origin at vortex center (m/sec) 

z Axial coordinate, along axis of vortex filament, with respect to local cylindrical 

coordinate system with origin at vortex center (m) 

5(x) Dirac (“delta”) function; x is the argument 

5'(x) Derivative of Dirac (“delta”) function, 5'(x) = d(5(x))/dx 

y Vortex filament circulation initial value, at t = 0 (m 2 /sec) 

v Kinematic viscosity (m 2 /sec) 

Q Rotor rotational speed (radians/sec) 

p Freestream density (kg/m 3 ) 

0 Polar coordinate, with respect to local cylindrical coordinate system, referenced to the 

vortex core center and the shear layer/slipstream boundary (radians) 

\j/ Stream function (m 2 /sec) 

\| /* Localized three-dimensional contribution to the vortex tangential velocity and 

vorticity, accounts for “line-sink” v z flow at the shear layer boundary (m 2 /sec) 

“Mean” stream function scaler value, function of r only, where V = (' ¥ + V* )d9 , 

(m 2 /sec) 

vpw Wake age (degrees) 

A\|/ w Net required wake age for vortex filament decay/collapse (degrees) 

co Flow vorticity 
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ROTOR VORTEX FILAMENTS: LIVING ON THE SLIPSTREAM’S EDGE 


Larry A. Young 
Ames Research Center 


SUMMARY 


The purpose of this paper is to gain a better understanding of rotor wake evolution in hover and 
axial flow by deriving an analytical solution for the time dependent behavior of vortex filament 
circulation and core size. This solution is applicable only for vortex filaments in the rotor far-wake. 

A primarily inviscid vortex/shear layer interaction (where the slipstream boundary is modeled as 
a shear layer) has been identified in this analytical treatment. This vortex/shear layer interaction 
results in decreasing vortex filament circulation and core size with time. The inviscid vortex/shear 
layer interaction is shown, in a first-order treatment, to be of greater magnitude than viscous diffu- 
sion effects. The rate of contraction, and ultimate collapse, of the vortex filament core is found to be 
directly proportional to the rotor inflow velocity. 

This new insight into vortex filament decay promises to help reconcile several disparate observations 
made in the literature and will, hopefully, promote new advances in theoretical modeling of rotor 
wakes. 


INTRODUCTION 


Even the most sophisticated of rotor free wake analyses neglect the decay of vortex filaments in their 
wake modeling. The assumption of constant vortex circulation and core size as adopted in free-wake 
analyses is too simplistic. And yet, to date, there is no new theoretical approach to replace this 
constant circulation and core size assumption. This dichotomy is noted in reference 1 . Most insights 
into vortex filament decay have their roots in fixed-wing aerodynamics. But, as the author will 
contend in this paper, application of these fixed-wing observations to rotorcraft is unjustified and 
ultimately leads to erroneous conclusions. 

What is not in dispute throughout the rotorcraft literature is that within three to five blade passages 
well-defined vortex filament structures can not be observed by means of any flow visualization 
technique. The analytical/computational treatment of the far-wake as evolving into an “ultimate” 
wake is, therefore, conceptually flawed. What has been the conventional wisdom, as regards this 
observed vortex filament decay in the far-wake, is that viscous diffusion effects are responsible for 
the process. Various turbulent vortex models, derived from fixed-wing analytical work, have been 
ad hoc modified and applied to attempt to account for, or match, the observed rotary-wing vortex 



filament decay process. There has not been a fundamental re-examination of the vortex 
dynamics/fluid mechanics for rotor trailed vortices. 

Other than qualitative flow visualization work, there is mixed empirical evidence for the time- 
dependent behavior of rotary-wing vortex filament circulation and core size. Most measurements 
have been made near, or before, the first blade passage. Correspondingly, these measurements have 
shown little wake age dependence of circulation and core size. Recent rotary-wing experimental 
observations, encompassing a greater range of wage ages than previously studied, are starting to 
appear in the literature. For example, one such study (A. Wadcock, Measurement of Vortex Strength 
and Core Diameter for a Model Tiltrotor, a soon to be published NASA CR) has presented hot-wire 
anemometry results that suggest that vortex circulation and core size both appear to decrease with 
wake age. This observation is counter to expectations based on the flow behavior of fixed-wing 
trailed vortices. However, an attempt to use the wide-field shadowgraph technique, reference 2, to 
measure rotor vortex core size as a function of wake age appeared to show the opposite trend from 
the Wadcock study. Reference 2 results implied that there was core size growth with wake age. 
Interestingly enough, both the Wadcock study and the reference 2 work conducted their vortex core 
size measurements with the same rotor/test-stand, in two different test entries. Both studies, though 
promising in their technique and initial results, should be considered as very preliminary bodies of 
work. 

Another noteworthy set of studies for vortex circulation and core size trends with wake age is the 
laser velocimetry work of references 3-5. However, though the authors of these works have made 
impressive progress in their measurement techniques, they have, to date, relied upon rather simplis- 
tic rotor systems that are not representative of actual aircraft. In particular, they obtained their test 
results for one-bladed, low tip speed rotors. As will be discussed later in the paper, the geometry and 
operating condition of rotor test models will have a substantial impact on the observed vortex circu- 
lation and core size trends. More systematic work using realistic scaled rotor models is required 
before empirical observations provide definitive answers to the rotary- wing vortex decay process. 

What is needed, more importantly though, is a fresh approach to the problem by reconsidering the 
fundamental differences between a trailed vortex filament from a fixed- and rotary-wing. The chief 
difference between the two types of trailed vortices, besides the helical geometry of the rotor vortex 
filament versus the rectilinear fixed-wing vortex, is that the rotor vortex is not surrounded by a 
uniform flow field as is the fixed-wing trailed vortex but is instead embedded in a shear layer. 

It is this flow feature, and the associated singular fluid mechanics, of a vortex in a shear layer that 
governs the decay process of a rotor trailed vortex filament. This area of investigation (vortex/shear 
layer interactions) has received very little analytical/computational work even though it would seem 
to be a fundamental type of fluid flow problem deserving considerable attention. 

This paper addresses the above noted oversight by developing an analytical model of vortex filament 
decay that accounts for vortex/shear layer interactions and is applicable to hover and axial-flight 
rotor wakes. 


This work is dedicated to the memory of Clara Jeannette McClelland Springer. 
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DESCRIPTION OF AN IDEALIZED FLOW PROBLEM AND 
ITS RELATIONSHIP TO ROTOR WAKES 


Figure 1 shows a vortex filament in the far wake region of a rotor slipstream. Sufficiently down- 
stream in the far- wake, the rotor slipstream flow is no longer accelerating and the boundary is no 
longer contracting. Further, the roll-up of vortex sheets inside the rotor wake has been completed 
and all vorticity is encompassed in the rotor trailed tip vortices. It is in the rotor far-wake that a 
number of simplifying assumptions can be made about the slipstream flow field, resulting in a 
“idealization” of the localized environment about a trailed vortex. First of all, the influence of wake 
contraction and nonuniform wake velocity profile will be neglected in this study. Further, the vortex 
filament helical angle and curvature can be assumed to be small and the flow problem can be treated 
as quasi-two-dimensional. Finally, the slipstream flow field is assumed uniform and the slipstream 
boundary is modeled as a discontinuous shear layer. 



Vortex Filament, 
refer to Figure 2a-c 
for cross-sectional 
close-up view 


v, 

1 


Rotor Disk 


Slipstream Boundary 


Figure 1. Vortex filament in the far-wake of a rotor slipstream. 


Figures 2(a-c) illustrate the vortex filament fluid mechanics (in a two-dimensional plane) at the 
slipstream boundary. Helmholtz’s Law requires that a vortex filament generated at the rotor disk 
accelerates to an equilibrium velocity of U + vi/2 along the slipstream boundary. The local cylin- 
drical coordinate system choosen for the idealized flow problem tracks along with an isolated vortex 
core center at this velocity of U + vj/2. Consequently, the velocity inside the shear layer is +Vi/2, 
while outside it the velocity is -Vj/2. This shear layer in the idealized flow problem is prescribed, 
time-invariant, and discontinuous. 
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Figure 2a - Radial velocity of unit 
volume of fluid just prior to 
crossing slipstream boundary, while 
within rotor wake. 




Figure 2b - Unit volume of fluid 
subjected to impulsive force when 
crossing slipstream boundary. 



Figure 2c - Radial velocity of unit 
volume of fluid just after crossing 
slipstream boundary; impulsive 
force results in momentum change. 

Figure 2. Problem description. 
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Figures 2(a-c) illustrate a “thought experiment” that asserts the existence of impulsive forces at the 
rotor slipstream boundary. Figure 2(a) depicts a unit volume of fluid prior to crossing the slipstream 
boundary, while it is still within the rotor wake proper; the unit volume has a radial velocity of +Vj/2. 
Figure 2(b) shows that the unit of volume of fluid, as it crosses the slipstream boundary, is subjected 
to an impulsive force. Figure 2(c) depicts the unit volume of fluid after crossing the slipstream 
boundary and just entering the freestream; the unit volume now has a radial velocity of-Vj/2. Con- 
sequently, the impulsive force applied to the unit volume of fluid as it crosses the slipstream has 
resulted in a net momentum change. This process occurs at 0 = 0, which is shown in figures 2(a-c), 
and is repeated at 0 = K (which is not shown). Force equilibrium is thus preserved as the impulsive 
force generated at 0 = 0 acts in opposition to the force at 0 = it. 

The vortex/slipstream shear layer interaction is assumed to be primarily an incompressible, inviscid 
flow phenomenon. It would seem, upon first consideration, that assuming inviscid flow while 
studying vortex filament decay does not make sense. Yet, it will ultimately be shown that the 
inviscid shear layer influence on vortex core size is much larger than viscous diffusion effects. 

Discontinuous shear layers, whether viscous or inviscid, are unstable (refs. 6 and 7). A rotor 
slipstream, or any naturally occurring shear layer, will consequently have a finite but steep shear 
flow gradient across its boundary. Nonetheless, a prescribed, time-invariant, discontinuous shear 
layer is very convenient as an idealization of a rotor slipstream. Because the shear layer is pre- 
scribed, and not part of the flow solution proper, the overall flow will not be subject to the 
instabilities noted in reference 7. 

It is the existence of these hypothesized impulsive forces, which are generated as vortical flow 
crosses the shear layer boundary, that makes rotor vortex filament decay uniquely and radically 
different from fixed-wing trailed vortex decay. The above thought experiment and corresponding 
idealization of the flow problem invoke fairly extreme simplifying flow assumptions. And, yet, the 
analytical solution of the described flow problem can and should be thought of as a “limiting case” 
for more physically realistic descriptions of the vortex/shear layer interaction. 

A Rankine vortex will be used to model the rotor tip vortex filaments in this paper. The analysis will 
assume that the Rankine vortex is isolated from the influence of other vortices along the slipstream 
boundary. A Rankine vortex was choosen for this analysis because of the simplicity of its modeling. 


SUMMARY OF THE APPLICABLE GOVERNING EQUATIONS 
FOR THE DESCRIBED PROBLEM 


Equations 1-4 are the governing equations for the vortex/shear interaction. Equation la is the 
incompressible r-momentum (cylindrical coordinate) Navier-Stokes equation. Equation lb is a 
derivative form of the r-momentum equation that has been integrated with respect to 0. Refer to 
appendix A for further details of the transforming equation la to lb. Equation lc explicitly includes 
the impulsive forces at 0 = 0 and 0 = n in the r-momentum equation. 


5 



d y r 

at 


+ v r 


d y r , v 9 ^ v r 
dr r 30 


+ v z 


dVj. 

dz 



r 


F + v 



2 3v e ^ 
r 2 dej 


I 


d_ 

at 


r 2n ^ 

J v r d e 

VO ) 



d y r , v 8 d y rl 
3r r 30 J 


d0 


» — 2v; 


3v 0 


dr 


10 = 0,71 


where F in equation la is given by 

F= -^f 


(la-b) 


(lc) 


The key modeling assumptions employed in equation lb are: the flow is inviscid; the flow is quasi- 
two-dimensional such that 3v r /3z = dv^/dz = 3v z /3z = 0, except that 3v z /3z = (l/v z )(3v z /3t) ^ 0 
when 0 = nrc, for n = 0,1,2,... ; the radial pressure gradient is primarily due to the vortex tangential 
flow pvg/r = 3P/3r for 0 * nit; the impulsive force can be related to the net momentum change 
across the discontinuous shear layer by 


-?- J ~{ f l[5(0) + 5(0 - 7t)] d 0} = -2vj 


3 vq 

3r 


l0=O,ir 


Equation 2a is an identity that follows straightforwardly from differential calculus. Definition of the 
two new parameters, vj/ and \j /*, eliminates 0 from the identity. Refer to appendix B for additional 
details as to the definition and application of equations 2a-b. 
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(2a) 


where 


\j) = J^ K (v|/ + \|/*)d0 (2b) 

Note that in equation 2b and in subsequent use in this paper, \|/* is the localized three-dimensional 
contribution to the vortex tangential velocity and vorticity, that accounts for nonzero v z flow at the 
shear layer boundary. 

Equation 3a is the general form of the vorticity conservation equation. Equation 3b is the derivative 
form of the conservation equation for a control volume about the vortex core boundary. This 
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integro-differential equation models the entrainment of negative vorticity into the vortex core; this 
negative vorticity is generated by the vortex/shear layer interactions. Refer to appendix B for details 
of transforming equation 3a to 3b. 
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The ¥* contributions to the vorticity conservation equation can be clearly seen in equations 3a-b. 
Note that V 2 (\j/ + V|/*), and not V 2 \j/, are used to define the vorticity. But, alternatively, 3\|//30 is 
used to define the vorticity convection velocity into the vortex control volume. 


Equation 4a is the incompressible z-momentum (cylindrical coordinate) Navier-Stokes equation. 
Equation 4b is a derivative form of the z-momentum equation. Refer to appendix C for the 
assumptions and approximations employed in reducing equation 4a to 4b. 
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(4a-b) 


Equations 1-4 are not dependent upon the type of finite-core vortex model employed in the analysis. 
The governing equations could be applied to any vortex model such as the Rankine or the Scully 
model (ref. 8). 


SUMMARY AND DISCUSSION OF A SOLUTION FOR A RANKINE VORTEX 
AS APPLIED TO THE DESCRIBED PROBLEM 


An analytical solution describing the decay of a Rankine vortex embedded in the shear layer of a 
steady, uniform slipstream is detailed in appendices A-C and summarized in equations 5-1 1 in the 
main body of this paper. 

The vortex/shear layer stream function (y + V) /*) that satisfies equations la-c has a spiral-like 
distribution in the region r c < r < r 0 (fig. 3). This spiral-like behavior of rotor vortex filaments is 
empirically well documented in the literature (e.g., refs. 9 and 10). Flow inside the rotor slipstream 
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Figure 3. Stream function (Y|/ + \|/*) (eqs. 5a-b). 


is shown in the lower half-plane of figure 3, and the flow outside the slipstream is in the upper half- 
plane. Equations 5a-b describe the vortex/shear layer stream function. 


V|/ = “VjSgn(9 - 7t)rsin0 + -^-f(t} 
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(5b) 
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(5c) 


Sgn(x) and Int(x) are the “signum” and “integer” functions, respectively. The constant k is 
proportional to the rate that the vortex core decays; the constant vj is the far-wake velocity, which 
is found from simple rotor momentum theory. Refer to appendix A for the details of the stream 
function derivation. 

The stream function, equations 5a-c, is valid for a finite radial region defined by 
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r c g(t)<r<— Vf(t) 

2trk 

The “finiteness” of the applicable flow region of the stream function is very important. It 
circumvents the prediction of physically unrealistic flow behavior far from the vortex core by 
confirming that the vortex/shear layer interaction is a “localized” phenomenon. 

The stream function, as written in equations 5a— c, is also subject to the flow constraint 


(^Slipstream Flow /^ r )/(^V’ Vortex /^ r ) > ® 

This flow constraint confirms that the stream function has been written to conform to the flow 
representative of rotary-wing slipstreams and vortices. Alternatively, bluff body wakes and shed 
vortices do not conform to this flow constraint. This issue will be discussed further later in the paper. 

A consequence of the derived stream function, equations 5a-b, is that negative vorticity is generated 
by the vortex/shear interaction (see eq. 6). 

co = k je - 7tlnt e + 71C + — 8'(0 - Jt) + — {28(0 - Jt)cos 0 + 5'(0 - Jt) sin 0} 


for r > r c g(t) (6) 

Equation 6 demonstrates that just as impulsive forces exist for a vortex embedded in a discontinuous 
shear layer, so does impulsive changes in vorticity. Most importantly, negative vorticity is created 
at the shear layer and will be shown to be convected/entrained into the vortex core. This result, 
equation 6, is crucial for the discussion that follows for the time-dependent behavior of the vortex 
filament circulation and core size — i.e., the decay process itself. Refer to appendix B for details to 
the vorticity expression, equation 6, derivation. 

The velocity and stream-function relationships used in this paper are subtly different from their 
conventional definition. They have been modified to account for quasi-two-dimensional flow. Quasi- 
two-dimensional flow has been defined in this paper as allowing localized, nonzero, time-dependent 
axial flow ( v z * 0 at 0 = 0 and 0 = Jt, whereas v z = 0 for 0 * n7t where n = 0,1,2,...) along the shear 
layer boundary. The velocity and stream function definitions have been modified such that now 
v 0 = +3(xj/ + v|/*)/3r and v r = -(l/r)(9\|//80). The tangential velocity expression has picked up an 
additional stream function term, \|/*, as compared to its conventional definition. 

The consequence of this quasi-two-dimensional flow revision to the stream function and velocity 
definition, when applied to the vortex/shear layer interaction problem, is that a constant radial flow 
into the vortex center is predicted, as well as an asymmetric tangential velocity distribution (see 
eqs. 7-8). The tangential velocity asymmetry can be seen in figure 4. The velocity profile shown 
in figure 4 is along the slipstream boundary (0 = 0 and 0 = Tt). Therefore, the asymmetry shown in 
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Tangential Velocity 



Radial Coordinate, r/rc 

Figure 4. Decay of vortices on a rotor slipstream boundary. 


figure 4 is distinct from the asymmetry encountered/expected when crossing the slipstream boundary 
(i.e., 0 = 7t/2 or 3 tc/2). 


v e = v i sgn(0 “ ft) sin 0 + ^ f ( t )|[“)[ 1 + s S n ( r " r c)]+ l 1 + s g n ( r c “ r )]J 
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v r = v i sgn(0 - 7t)cos0 — k 


( 8 ) 


Derivation of radial and tangential velocity analytical expressions follow directly from the derived 
stream function, equations 5a-b. 

Radial flow and vorticity entrainment into the vortex core causes the vortex circulation to be time 
dependent. Vortex circulation is found to be a decreasing function of time and wake age. Refer to 
equations 9a-b and figure 5. 
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Circulation = yf (t) 


(9a) 


where 



time, t 

Figure 5. Vortex circulation and core size as a function of time (eqs. 9 and 10). 


y and r c are the initial values, respectively, of vortex circulation and vortex core radius. The function 
g(t) is defined in equation 10b. The function f(t) encompasses the time dependency of vortex circula- 
tion. The function g(t) reflects the time dependency of the vortex core size. Refer to appendix B 
for the details of the derivation of the above analytical expression for the vortex circulation decay. 
The rate of decay of circulation can be seen to be directly proportional to rotor inflow velocity. 
Identification of a potential flow mechanism that results in decreasing circulation is quite surprising. 

Even more surprising is the analytical result for the time dependence of the vortex core size. Core 
size is also found to decrease with time and wake age. Refer to equations lOa-b and figure 5. 

Core Radius = r c g(t) (10a) 


where 


yf(t) + tt 2 kr c g(t) = Jt 2 kr c + y - (2vjk + tt 2 k 2 Jt (10b) 

rc is the initial value for the vortex core radius; equation 5b provides the value of the constant, k. 

Equations 9b and 10b are a nonlinear, simultaneous set of equations that are iteratively solved for 
f(t) and g(t) for each time, t. Refer to appendix B for the details of the derivation of the solution of 
the vortex circulation and core size time dependent behavior. 
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The prediction of vortex core size contraction is contrary to conventional thinking, which is based in 
large part upon observations made from rotary-wing flow visualization studies. This discrepancy 
will be discussed further later in this paper. But a key consideration in interpretation of rotor wake 
flow visualization results is the assumption of the relative magnitude of axial flow, v z , along the 
vortex filament length. 

Work from this paper (eqs. 1 la-b) predicts large vortex filament axial velocities along the shear 
layer boundary. 


v = e ^k8(0-nn)t/ r 


for t«A\|r w /Q (11a) 


v z — > Ae sgn ( 0-7 ^ rsin0 for t->~ 


(lib) 


where 


A * 



+ v 2 j-2In(v) 


(11c) 


Refer to appendix C for details of the axial flow solutions. Equation 1 la describes a “line-sink.” A 
portion of the mass flow crossing the slipstream boundary is redirected and convected axially along 
the shear layer. Equation 1 la is a limiting case for more realistic shear layer models (where shear 
layer gradients are steep but not discontinuous), and so equation 11a suggests that axial velocities for 
rotor tip vortices are large, but finite, and increase with time. Ultimately, at t — > <», the axial velocity 
distribution of the asymptotic solution, equation 1 lb, suggests that a swirl distribution develops at 
the slipstream boundary that counterbalances the swirl in the rotor slipstream from the blade root 
vortex. Through equivalence of kinetic energy (an energy balance between the swirl from the blade- 
root vortex versus that of the swirl induced from the trailed tip vortices), the constant A was defined, 
equation 11c. 


RE-EXAMINATION OF THE ANALYSIS ASSUMPTIONS IN LIGHT 
OF THE DERIVED SOLUTION 


It might be supposed that reliance on an inviscid analysis could result in invalid conclusions as to 
vortex filament circulation and core size time-dependent behavior. A first-order examination of the 
relative contraction and expansion rates of inviscid versus viscous effects will show that the inviscid 
vortex/shear layer interaction dominates the fluid dynamics of rotor vortex filaments. Examining 
more closely this assertion, i.e., 
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For t » 0, the predicted inviscid vortex core contraction rate will be compared to the viscous 
diffusion rate of a Lamb-Oseen laminar trailed vortex (ref. 1 1). Or 
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Using typical values for vj and v for equations 12a-c, the inviscid contraction rate is estimated as 
being a couple of orders of magnitude greater than the laminar vortex core expansion rate. For 
example, for A\(/ w > 2 degrees (t » 0), the laminar expansion rate is less than 0. 1 foot per second. 
But, for a hovering rotor, for 0.01 < Ct < 0.015, which is typical of 1G thrust loaded rotors, the 
predicted inviscid contraction rate is 50 to 60 feet per second. Some authors have suggested that the 
contribution of turbulent viscosity to a trailed vortex diffusion rate may be as much as an order of 
magnitude greater than the laminar estimate (see ref. 3). Nonetheless, if this supposition was proven 
correct for rotor trailed vortices, the inviscid vortex/shear layer interaction contraction rate is still 
greater than this inflated viscous diffusion rate estimate. Therefore, the inviscid flow assumption 
employed in the beginning of this paper appears to be, ultimately, justifiable. 

Another assumption employed in this paper’s analysis that bears further scrutiny is that the radial 
pressure gradient is dominated by quasi-steady vortical flow. That is, it has been assumed that 

= (13 

r p dr 


The static pressure can be separated into two components: a quasi-steady pressure component that 
is dominated by the vortex flow and an unsteady, non-vortex-driven pressure component from the 
vortex/shear layer interaction. That is: P = Pqsv + PuNV- The implication of this observation is that 
the above assumption (eq. 13) is equivalent to the assertion represented by equations 14a-c (for 
0 * 0 * tc where the pressure gradient contribution to the surface forces has been encapsulated in 
the impulsive shear forces, 1(5(0) + 5(0 - tc)), which have been explicitly separated out): 
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Therefore, if equation 13 is true then, given equations 14a-c, the inequality constraint can be 
re-written as 
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Or, making substitutions from equation 5a-b, 

+ kje - K Int^— » 0 


(15a) 


(15b) 


Using typical values for f(t), k, r c , and y the left hand side of the equation 15b inequality is found to 
be quite large — one to two orders of magnitude greater than unity. Using the reference 2 values for 
f(t), k, r c , and y the ratio r/r c has to be greater than the value of 4.2 before the magnitude of the left 
hand side of the inequality drops below unity. The stream function is valid to only ro/r c = 4.3. The 
flow assumption represented by equation 1 3 thus appears valid. 

Besides viscous flow effects, the distributed vorticity of real rotor vortices versus the concentrated 
vorticity in the core of a Rankine vortex and the interactional influence of multiple rotor tip vortices 
on each other will each tend to counterbalance the predicted core collapse. Additionally, the 
influence of rotor slipstream acceleration on the vortex circulation and core size time dependent 
behavior needs to be assessed in further work. 


IMPLICATIONS OF THE VORTEX/SHEAR LAYER INTERACTION SOLUTION 
WITH RESPECT TO OBSERVATIONS MADE IN THE LITERATURE 


The predicted flow characteristics from this paper’s work provide insight into several puzzling and 
conflicting observations made in the literature regarding rotor wake behavior. Using the analytical 
results from this paper, many of the these observations can at least be partially reconciled with each 
other. 

As noted earlier, a key empirical observation for rotor wakes, as reported innumerable times in the 
literature, is that the rotor far wake loses, after three to five blade passages, its well-defined vortex 
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filament structure. Previously, the decay mechanism for the far-wake trailed vortex filaments has 
been assumed to be due to viscous diffusion of the filament vorticity. This paper has promoted an 
alternate hypothesis and intepretation of the filament decay mechanism and, further, has derived a 
solution for this inviscid vortex/shear layer interaction that predicts vortex circulation and core size 
to decease with time/wake age. 

But, as noted earlier, the prediction of vortex core size contraction is contrary to conventional 
thinking. Most experimental studies reported in the literature on rotor wakes and vortex core size 
measurements have relied on smoke or particle flow visualization or, alternatively, velocity 
measurements near the rotor disk plane. The vortex core has been equated, or assumed directly 
proportional, to observed “particle voids” in vortex filaments for seeded rotor wakes. And, indeed, 
these vortex filament particle voids, for rotor wakes, do grow in size with time/wake age. However, 
the author asserts that smoke flow visualization is an invalid measure of vortex core size. 

The fascinating flow visualization photographs of reference 12 are the clearest illustration that the 
smoke particle void is not the same as the vortex core. In the reference 12 photographs, smoke flow 
visualization results (i.e., the particle void) are observed simultaneously with fortuitous images of 
the vortex core as seen by means of condensation. In many cases, the smoke particle void is much 
greater in size than the vortex core. 

Smoke flow visualization results are an invalid measure of vortex core size because of the influence 
of large axial velocities, v z , existing along the vortex filament length that sweep regions greater than 
the vortex core clear of smoke particles. There is considerable analytical and experimental work that 
confirms the existence of large axial velocities along a trailed vortex filament length. The prediction 
of nonzero axial velocity along the vortex filament length, as derived in this paper, qualitatively 
agrees with helium bubble flow visualization observations recounted in references 1 and 13. Finally, 
the rotor swirl at the slipstream boundary that is implicit in the asymptotic axial velocity solution, 
equation 1 lb, matches qualitatively rotor wake swirl measurements noted in reference 10. Refer to 
figure 6. 



Nondimensional Blade 
Radial Station 


Figure 6. Swirl distribution of the rotor far wake, accounting for the asymptotic filament axial flow 
resulting from the inviscid vortex/shear layer interaction. 
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On the basis of both measurements and predictions, axial velocities are clearly of a magnitude to 
invalidate the assumption that the particle void at the center of a vortex is representative of the 
vortex core radius. The particle void as seen in smoke visualization work should asymptotically 
approach the value of r o (0) — where r 0 (t) = yf(t)/27t 2 k and so r o (0) = y/2^ 2 k — rather than r c g(t). 
This hypothesis can be tested in future studies by comparing predicted r o (0) values to particle void 
trends. (Reference 13 is one of few studies in the literature to provide quantitative data of the particle 
void size trends with wake age. Unfortunately, reference 13 did not examine large enough wake ages 
for an asymptotic limit to be observed.) If the above hypothesis is true, i.e., the particle void is 
proportional to the r o (0), the asymptotic wake age limit for particle void size will be dependent upon 
tip speed and blade count for hovering rotors. 

The apparent discrepancies in vortex core size trends between the work of A. Wadcock (Measure- 
ment of Vortex Strength and Core Diameter for a Model Tiltrotor, a soon to be published NASA 
CR) and the shadowgraph results of reference 2 can be commented on in light of the predictions 
of this paper. The study by Wadcock made hot-wire anemometry measurements of vq and was, 
therefore, able to directly estimate core size. The shadowgraph results of the reference 2, on the 
other hand, would tend to be biased towards a measure of the resultant velocity vector, i.e., 

•\/ v z + v e + v r * an d not just V0. Interestingly enough, smoke flow visualization and shadowgraph 
test results will both be adversely impacted as a measure of vortex core size by not being able to 
account for non-negligible axial velocities, v z , along the filament lengths. 

The majority of research into rotor vortex circulation and core size trends with wake age have relied 
upon small-scale, low thrust, and low tip speed rotors. Also, one-bladed rotors have often been tested 
to keep the wake geometry as simple as possible. See references 3 and 13, for example. There is a 
critical scaling issue, though, in the use of small-scale, low tip-speed rotors and one-bladed rotors for 
studying vortex decay mechanisms and generalizing the results to rotors more representative of 
flight hardware (refer to eq. 16). 
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(16) 


where A\|/ w is the net amount of wake age required for decay of rotary-wing vortex filaments. Refer 
to appendix B for details of the derivation of equation 16. The derivation of this wake age expression 
for vortex decay follows directly from the nonlinear expressions for the vortex core size and 
circulation, equations 9a-b and lOa-b. 

Examining equation 16 closely, the equation suggests that one-bladed rotor test results can 
potentially provide misleading insights into vortex circulation and core size trends. The vortex 
filament decay of the rotor studied by A. Wadcock (Measurement of Vortex Strength and Core 
Diameter for a Model Tiltrotor, a soon to be published NASA CR) — a small-scale rotor but with 
typical blade count and tip speeds for full-scale rotorcraft — is predicted to occur over roughly 
A\|/ w = 60-80 degrees. The reference 3 rotor (low tip speed and a single rotor blade) theoretically 
takes five times longer to decay out, i.e., A\j/ W = 300-400 degrees. Reference 3 identified a gradual 
decrease in circulation with wake ages in the far wake region; mixed results were seen for the core 
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size trend — though the authors ultimately concluded that there was overall growth in core size with 
time. Because of the intrinsically lower vortex decay rate for the reference 3 rotor, these results can 
not be considered conclusive. 

Figures 7 and 8 provide a more detailed examination of the relative influence of blade count, thrust 
loading, and axial-flight (or climb) advance ratio on the predicted decay rate, which is inversely 
proportional to A \]f w , for the vortex core collapse. Both figures are based on predictions from 
equation 16. Figure 7 shows the influence of thrust and axial-flight advance ratio on A\j/ W . Both 
rotor thrust and axial-flight advance ratio have a profound effect on the decay rate of trailed vortex 
filaments. Increasing axial-flight (or climb) advance ratio decreases the decay rate. An interesting 
interpretation of figure 7 is that as the axial-flight advance ratio asymptotically approaches infinity 
the rotary-wing vortex filaments will appear and behave almost like fixed-wing trailed vortices. At 
high axial-flight advance ratios the inviscid vortex/shear layer interaction becomes negligible and 
viscous diffusion will dominate the vortex decay process just as it does for fixed-wing trailed 
vortices. Figure 8 shows the influence of blade count and thrust for hovering rotors on A\j/ W . In 
hover, rotor thrust has little impact on vortex decay rate, but increasing blade count and tip speed 
substantial increases the decay rate. 

An interesting question is what happens to the vortex core size and circulation trends if the vortex 
tangential flow is counter to the slipstream flow, i.e., (^slipstream Flow/^ r )/(^¥vortex/^ r ) < 0. 

This type of vortex/shear layer interaction is now representative not of the flow field of a rotor wake 
but is instead applicable to bluff body/shed vortex flow. It is asserted that equations 5-11 will also 
hold for bluff body /shed vortex flow — if the analysis is implemented with a change in sign of vortex 
circulation and the constant k (but not the value of the constant). 



V T 


Figure 7. Vortex decay rate as a function of thrust and axial flow advance ratio (three bladed rotor). 
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N 

Figure 8. Vortex decay rate as a function of thrust and blade count for a hovering rotor. 


To test this assertion, the results of this revised implementation of equations 5-1 1 will be qualita- 
tively compared to the results of a computational fluid dynamics (CFD) treatment of the bluff body 
vortex/shear layer problem. Reference 14 is one of the few analytical or computational treatments 
of the vortex/shear layer problem. In the particular case of reference 14, a vortex embedded in a 
viscous, constant shear flow gradient was studied with an unsteady, Navier-Stokes CFD code. Given 
the reference 14 definition of the vortex/shear layer Reynolds number, the discontinuous, inviscid 
shear layer studied in this paper would yield Re — » A qualitative comparison can be made 
between this paper’s analysis and the reference 14 results by studying the reference 14 high 
Reynolds number cases (Re > 1000). In both cases, vortex circulation increases rapidly in strength 
with time. No predictions were made in reference 14 for core size trends; equations 5-11, however, 
predict a slow growth of core size with time. 

Despite this initial attempt to reconcile the experimental work of the past with the results of 
this paper, additional experimental investigations are necessary to corroborate and extend the 
analytical work described in this paper. Work in this area, concentrating on nonintrusive flow 
measurements that survey a large range of wake ages, is critical to understanding vortex filament 
decay mechanisms for rotors. Acquiring test data for more realistic rotor models is also crucial. 


CONCLUSION 


Many rotary-wing analyses rely upon the assumption of constant vortex circulation and core size in 
their formulation. However, an analytical treatment of vortex/shear layer interaction suggests that 
there is a substantial time (or wake age) dependence of vortex circulation and core size. Both are 
shown to decrease with time. The results of this paper will, hopefully, promote new insights into 
rotor wake flow behavior. 
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APPENDIX A - STREAM FUNCTION 


Derivation of the Governing Equation 


Beginning with the three-dimensional, incompressible r-momentum Navier-Stokes equation for a 
vortex embedded in a shear layer (cylindrical coordinates, with the origin centered on the vortex 
filament axis, and z as the coordinate along the vortex filament axis; the shear layer boundary 
assumed to be locally straight and aligned along 0 = 0,7t): 


^L +v ^r+Ie^L + v 
9t r 3r r 30 z 3z 


^L-^ = F + V [V 2 V r 


V r 2 ^ 

r 2 r 2 30 j 


(Ala) 


F = --? + gr + -|-{i[8(e) + 8(e - *)]} (A lb) 

p 3r p dr 1 

An impulsive force term has been incorporated into the r-momentum Navier-Stokes equation, 
equations Ala-b, as a “driving” force that results from vortex flow crossing the slipstream 
boundary. 1(8(0) + 5(0 - 7t)) is an impulsive force function where 5(0) and 5(0 - n) are Dirac 
functions (e.g., ref. 15). (In short, 8(x - x 0 ) = 0 for x * x 0 and 5(x - x G ) = °° for x = x G ; finally, 

J°° 8(x - x D )dx = 1 .) The impulsive force accounts for the net radial momentum change that occurs 
asvortex-driven flow crosses a prescribed, time-invariant, and discontinuous shear layer. (A “free” 
shear layer, on the other hand, would not remain discontinuous for large periods of time under the 
influence of viscosity (see ref. 6).) Refer to figures 2a-c. An expression for the vortex/shear layer 
impulsive force function, 1(5(0) + 5(0 - ji)), will be defined later in this appendix. 


It is assumed that the decay of a vortex filament embedded in a shear layer is primarily a potential 
flow phenomenon. Additionally, the flow is assumed quasi-two-dimensional such that 
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Therefore, invoking the inviscid, quasi-two-dimensional flow assumptions reduces equations Ala-b 
to 


av L + v 5v r [ v e 5v r v e 
5t r 5r r 50 r 


F 


(A2) 


Finally, it is assumed that, for vortex dominated flows, the following approximation can be made: 

vq __ 1 5P 
r p 5r 


The above is only strictly true for steady vortex flow fields and, therefore, inherent in this 
approximation is the assumption that the vortex core contraction/expansion rate is small in 
magnitude. Proceeding with this approximation 


d v r y d y r v 6 ^ v r 
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Integration of equation A3 with respect to the angular coordinate, 0, transforms the impulsive force 
to an incremental momentum change across the slipstream boundary (see ref. 15). This integral 
operation against equation A3 is performed, drawing upon an analogy from applied mechanics as to 
the mathematical treatment of impulsive forces in classic dynamic system analysis. Integration of an 
impulsive function over the complete domain of the function’s independent variable yields a finite, 
net momentum change. In classic dynamic problems, the impulse function independent variable is 
time and the function integrated across zero to infinity yields the dynamic system’s net momentum 
change. In the case of a vortex interaction with the slipstream shear layer, the impulse function’s 
independent variable is 0 and so integrating from 0 < 0 < 27t yields the total net momentum change 
for a unit volume of vortical flow as it crosses the slipstream boundary twice in one revolution. 
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Or, making the substitution for the incremental momentum change, i.e.. 
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where the term dv e /dr| e=0 denotes the derivative of the tangential velocity with respect to the 
radial coordinate as evaluated at either 0 = 0 or n, as the derivatives at these two polar coordinates 
are equal in magnitude and sign. Therefore, in interpreting equation A5, it should be noted that as a 
unit volume of vortical flow crosses the slipstream boundary the flow sees an impulsive step in 
radial velocity of v;, towards the core center. 


Solution for a Rankine Vortex 

In order to proceed with the analysis it is now necessary to assume a general expression for the 
stream function that is based, in part, on unknown functions that will later be solved for. It is 
assumed that there are three major components to the stream function for the flow region defined by 
0 * nn, where n = 0,1,2,.... A separate stream function term, y*, will defined later for 0 = rm to 
completely describe the quasi-two-dimensional flow regime. Letting 

V - ¥ Uniform Slipstream Flowfield + V Rankine Vortex Flowfield + ¥ Vortex/Shear Layer Interaction 


where 


¥ Uniform Slipstream Flowfield 
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v; sgn(0 — 7t)rsin0 


Note that if the above expression is transformed from polar coordinates to Cartesian coordinates, and 
the appropriate partial derivatives are taken to define the flowfield velocities (v x and Vy versus v r 
and vq), it can be shown, as expected, that 
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and, further. 


¥ Rankine Vortex Flowfield ^ ^(0' 
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where, for t = 0, by definition, the unknown functions have the value g(0) = f(0) = 1 . Therefore, for 
t = 0, the above expression reverts to the stream function — and correspondingly the velocities (v r 
and vq) — of the classic analytical expressions for a Rankine vortex, i.e.. 
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which is comprised of the product of three unknown (but assumed separable) spatial and time- 
dependent functions. 

Therefore, in summary, it is assumed that the stream function can be given by the expression 
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+ h(t)a(r)x(9) 


where 


Circulation = yf(t) (A7a) 

Core Radius = r c g(t) (A7b) 

Equations A6 and A7a-b represent the superposition of a shear layer and a Rankine vortex (modified 
to incorporate time dependent functions for circulation and core radius) and the addition of a vortex/ 
shear layer interaction term. See reference 16, for example, for additional discussion as to Rankine 
vortices. A Rankine vortex was selected to model this problem as it is a finite-core vortex model that 
inherently satisfies the steady-state potential flow equation. It is recognized that a Rankine vortex 
has substantial limitations as to its accuracy in modeling the velocity distributions of rotor trailed 
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tip vortices. Yet, for the purposes of this paper, which emphasizes deducing the first-order, time- 
dependent behavior of vortex core size and circulation, this should not pose a concern. 


To accommodate the assumption of quasi-two-dimensional flow, the definition of velocity terms 
with respect to the stream function will be subtly modified: 

v '=4f bM ves l (M ' +v * ) 

(A8a) 

where 


M i* = e(t)<p(r)p(0) 

(A8b) 

p( 9 )- m „{T) 
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Note that p(0) is a variation on the step (“Heaviside”) function— except that the step is progressive 
with each crossing of the slipstream/shear layer boundary. p(9) is a generalized function just as the 
step, signum, and the Dirac (“delta”) functions are (ref. 17). In accordance with the definition/ 
assumption of quasi-two-dimensional flow, it can be seen that ^ v 9/^®|0 =nlt 06 0(p(®))/00 where, as 
earlier hypothesized, 3(p(0))/90 «= 8(0 — nrt). The coefficient b in equation A8c is a constant. 

Therefore, the velocity functions that satisfy equation A6 must have the form: 

v r = -jvj sgn(0 - 7t)cos0 - -j-h(t)a(r)-^(x(0)) (A9a) 
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2 471 r c g(t) j (A9b) 

+ h(t)x(0) (a(r)) + e(t)p(6) j- (tp(r )) 

dr dr 

where sgn(x) is the signum function, and sgn(x) = -1 for x < 0, sgn(x) = 0 for x = 0, and 
sgn(x) = +1 for x > 0. 

The signum and integer functions can be expressed in terms of unit step or Heaviside functions. 

Dirac (delta) functions, 8(x), have already been employed in the derivation of the governing equa- 
tion, equation A5; Dirac functions will also manifest themselves in the evaluation of equation A5. 
The step and Dirac (delta) functions are “generalized” functions (see ref. 17). Formulas for integral 
and differential operations applied to generalized functions are well established. 
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The identification of analytical expressions for f(t), g(t), h(t), a(r), x(0), e(t), and cp(r) that satisfy 
equation A5 (shown immediately below, once again, for convenience) is the principal undertaking of 
this paper. 
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Substitution of equations A9a-b (the assumed functions for the radial and tangential velocities) into 
equation A5 yields the following results for the integro-differential terms in equation A5: 
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Flow visualization of rotor wakes — showing the interaction of the tip vortices with the slipstream — 
reveals spiral-like flow patterns (e.g., refs. 9 and 10). This leads to the insight that the vortex/ 
slipstream interaction term in the stream function could perhaps be approximately represented by 
the simple analytic expression for a spiral, t|/ vortex/Shear Layer Interaction oc ^* 0 * 

Proceeding with this insight/hypothesis, it will be shown that equation A5 (for r > r c g(t)) can be 
satisfied if the following functional forms for a(r), y(9), and (p(r) are employed: 


And, further, assume that 


a(r) = ar (Alla) 

*(0) = c0 (Allb) 

cp(r) = a(r) = ar (Allc) 

e(t) = h(t) (Alld) 
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The assumptions explicit in equations A1 lc-d can be rationalized on the basis of functional analysis. 
Substituting the assumed functions, equations A1 la-d, into equation AlOa-d, for a(r), %(©), £(t), 
and <p(r), yields: 
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Equations A12a-d can be further reduced by noting the following identities (ref. 17): 


— (sgn(0-7t)) = 28(0-7t) 


(A13a) 
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This is the well-known “sampling property” of the Dirac (delta) function, /(x) is any continuous 
function. 


EL, /( x ) u ( x - x o) 5 ( x - x o) dx = J_^/( x ) k5 ( x - x o) dx = k /( x o) (A13c) 

Where u(x - x 0 ), as used in the above equation, is the unit step function, /(x) is any continuous 
function, and k is a constant (where 0 < k < 1). Various authors assign several different values to k 
but the value used by the majority, and adopted in this paper, is k = 1/2 . 
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(A 13d) 
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Applying equations A13a-g to equations A12a-d gives 
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The constant b, in the function p(0), can be determined by taking into account symmetry 
considerations for the vortex interaction, i.e.. 


v ^le=7t + " v ole=o + 


which gives (making the appropriate substitutions from eq. A9a into eq. A15) 
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Substituting the integro-differential terms of equations A14a-d back into equation A5 gives 
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Equation A17 can only hold true for r > r c g(t). Therefore, 
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Equation A18 then requires that 


dt 


(h(t)) = 0 


or, more specifically, 
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Equation A19 will only be valid until such time as the vortex completely decays out. Once the 
vortex decays out then h(t) = 0. 


To recap, the stream function that satisfies the r-momentum equation for a Rankine vortex embedded 
on a slipstream shear layer has been found to be (for r > r c g(t)) 
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where k, a constant, is given by 
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The positive sign of the constant, k, is consistent with the inward (towards the vortex core) radial 
velocity that the impulsive force imparts across the shear layer. Given equations A9a-b, A 16, and 
A20a-b, the velocities vq and v r are therefore 


v r = -^ Vj sgn(0 - 7t)cos0 - k 


(A21a) 


v 0 = - ^ v i sgn(0 - 7t) sin 0 + ^ f (t)|(j |l + sgn(r - r c )] + l 1 + s § n ( r c “ r )] j 




+ k 


0 - Ttlnt 


V 


m 


(A21b) 


Unfortunately, a solution for the stream function for the flow in the vortex core does not suggest 
itself at this time. This limitation in the analytical treatment of the vortex/slipstream problem does 
not adversely impact the evaluation of the vortex core size and circulation time dependence, as can 
be seen in appendix B. 

Finally, it should be noted that there is an outer radial limit, r Q , for which the derived stream function 
applies. Beyond this limit there is no longer any mass transfer (via nonzero vortex tangential 
velocities) across the shear layer. This outer limit can be derived by noting that 

v ele=o = v ele=K = 0 at r = r ° 


which, given equation A21b, yields 


— -^^-7tk = 0 
2k r 0 


Therefore 


r 0 =— T-f(t) (A22) 

2n k 

In conclusion, equations A20-A21, inclusive, are valid for the region defined by the radial limits 

r c g(t)<r<— ^-f(t) 

2f 7U 1c 
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APPENDIX B - VORTEX CORE SIZE AND CIRCULATION TIME DEPENDENCE 


Next, this paper will concentrate on deriving expressions for f(t) and g(t) for the stream function for 
r > r c . It is possible to deduce the functional dependence of the vortex core size with time by noting 
the following relation. 


dr 

3t 


1 


cKj/ 
cKfO 3t 
9r 


(Bla) 


where 


V = Jq (V + ¥*) d0 


(Bib) 


or, making the appropriate substitutions, 




( „ V 




— 7t 2 kr 


Refer to appendix A for the definitions of V| / and \|/*. The expression for \j> requires that 
equation Bla reduce to 




f(t) 27t 2 r c k 
g(t) Y 


at gW 


Or, alternatively, using the chain rule for derivatives, 


— -41 
dg 


f 


u 


2n v r k 


C*v _ 


= 0 


(B2) 


The above is a homogeneous, linear ordinary differential equation which has the solution 


f(t) = 


^2Tt 2 r c k^ 
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g(t) + 


/ " ] _2jt 2 r c k^ 


V 


3y 


g(tr 


(B3) 


) 


The initial value condition (g(0) = f(0) = 1) has been accounted for in equation B3. 
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To arrive at final closure for the functions f(t) and g(t), it is necessary to evaluate the vorticity 
conservation equation (ref. 18). 


^-(jcodvj = -|(con)udS + J(n xF)dS (B4a) 


where 


co = -V 2 vj/ 


(B4b) 


and 


V 2 = 



i _a i a 2 

r dr r 2 30 2 


Note that V 2 (\|/ + \j/*) and not V 2 v|/ is used to defined the vorticity, but 3y/d9 is used to define 
the vorticity convection velocity into the vortex control volume instead of 8(\|/ + \|/*)/80. Refer to 
appendix A, equations A20a-c, for the derivation of and expressions for \| / and vj/* . The vorticity 
conservation equation, for a control volume about the vortex core, is reducible to (noting that there 
are no normal surface forces at the vortex core boundary): 


d_ 

8t 



2n 

I 


j^ + r 

[39 c 


2 d 


~ g(t) }■ * ( V + V* )d0| 


(B5) 


The vorticity convection into the vortex core, due to the shear-layer-induced radial flow, has to take 
into account the vortex core contraction/expansion rate. The vortex core contraction/expansion rate 

expression is r c — g(t) and has been included in equation B5. 
dt 

Given equations A20a-c and A21a-b, the vorticity can be defined by noting that 

V 2 (\j/ + i|/*) = -kje - - ^6'(0 - n) - ^-{28(0 - 7t)cos0 + 8'(0 - 7t)sin 0} 


for r > r c g(t) 


(B6) 


In the above nomenclature 8'(x) = d(8(x))/dx. Making the substitution of equation B6 into B5 and 
noting the following three points: (1) the control volume boundary is a function of time; (2) only the 
vorticity generated by the vortex/shear layer interaction for r > r c g(t) is convected into the control 

~ K ( ) 

volume; and (3) i r u V 2 \|/rdrd0 = yf (t) for the vorticity contained in the vortex core. Therefore 
0 
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(B7) 


y-j-f(t) = k J je- rclnt ^^— j - tc8'(0 - 7 t)||-^ VjSgn(e - 7t)cos0 + k + r c •^•g(t)|d0 

-Vj J |[28(0-7c)cos0 + 5 / (0-7i)sin0]}|-^-v i sgn(0-7c)cos0 + k + r c ^g(t)|d0 

The following integral relations are applicable in evaluating equation B7 

f /(x)8'(x-x 0 )dx = -/'(x 0 ) (B8a) 

J — oc 


/(x) u(x - x 0 )8'(x - x 0 )dx = J_^/(x)k5'(x-x 0 )dx = -k/'(x 0 ) (B8b) 


where u(x - x Q ), as used in the above equation, is the unit step function; /'(x) is the derivative of a 
continuous function, /(x), which in itself is also continuous; and k is a constant (where k is arbitrary 
in value). 


Therefore, application of equations B8a-b to equation B7 results in a first-order, linear differential 
equation with the two unknown functions, f(t) and g(t) 


Y^-f(t) + (2vj +tt 2 r c ^g(t)lk + 7t 2 k 2 =0 
dt V dt J 


(B9) 


Integrating equation B9 with respect to time yields 

yf(t) + 7t 2 kr c g(t) = 7t 2 kr c + y -^2vjk-i-7t 2 k 2 jt 


(BIO) 


An initial value constraint of g(0) = f(0) = 1 has been accounted for in equation B 10 by solving for 
the integration constant. The above equation is a nonlinear implicit solution for the vortex core size. 
The vortex core radius can readily be seen to decrease with time by noting an approximation for g(t) 
at very large units of time when f(t) and, therefore, the circulation approaches zero. 


g(t)«l + 


Y 

7t 2 kr c 


/ 2v i k + Tt 2 k 2 ^ 
v rt 2 kr c 


(Blla) 


as 


7t 2 kr c + y 
2v;k + 7t 2 k 2 


(Bllb) 
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The interdependence of the initial value of circulation and the slipstream inflow velocity for rotor 
wakes can found by simple momentum/vortex theory for rotors. The trailed circulation for a tip 
vortex filament shed in rotor axial flow can be found by the expression (see ref. 8): 


Y = 


4ttRv 

NV t 


(U + v) 


(B 12) 


where N is the number of blades and Vj is the rotor tip speed. The mean inflow at the rotor disk, v, 
is 


v — 



+ 


— V- 
2 


(B13) 


As this analysis is applicable for the mid- to far-wake for rotors (i.e., the flow is steady, 
nonaccelerating), then based on simple momentum theory for rotors 

V; = 2v (B14) 

To conclude this discussion an expression for the wake age required to completely collapse the 
vortex will be derived. Referring to equation B1 lb 


A\j/ W = £2t = Q 


rc 2 kr c + y 
2v ; k + 7i 2 k 2 


(B 15) 


Making appropriate substitutions from equations B12, B13, and B14 into equation B15 yields 


Ay w = 


^ n 2 + 8 ^ 




V v ) 


r 2rc 2 + 16 ^ 
TtN 


i + H 

V 


(B16) 
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APPENDIX C - VORTEX FILAMENT AXIAL FLOW 


Axial Flow Along Slipstream Boundary for t « A\|/ w /Q 

The vortex/shear layer interaction problem has been assumed to be quasi-two-dimensional. This has 
been defined to be 


3v 

— — = 0 for 9 n7t where n = 0,1,2,... 

dz 


but 


dv z _ dv z dt 


1 dv 


dz 


3t dz v z 9t 


— * 0 for 9 = nrt 


The continuity equation along the shear layer/slipstream boundary can be rewritten as 


1 3 y z _ [ 9v r f v r , 1 dv 9 
3r r r 39 


v z 3t 


which gives 


v 


2 





(Cl) 


(C2) 


Finally, expressions for v r and vq have been derived in appendix A and are given by equations 9 and 
10 in the main body of the paper. Making the appropriate substitutions 

v z =e” k5(0 ' n,t)t/r (C3) 

Equation C3 predicts that v z is a constant for r < r 0 (refer to appendix A, eq. A22, for the 
definition/derivation of r Q ) for all 9 — except along the slipstream boundary, i.e., 9 = nTt. For 
9 = mt, equation C3 describes a line-sink; some of the tangential mass flow crossing the slipstream 
boundary is redirected and convected axially along the shear layer. Continuity has, therefore, 
been shown to be preserved even at the discontinuous slipstream boundary (if one overcomes the 
conceptual hurdle of an infinite axial velocity, acting along an infinitesimal shear layer thickness, 
that increases with time). Equation C3 should be seen as a limiting case for more realistic shear layer 
models (shear layer gradients that are steep but not discontinuous). Therefore, axial velocities for 
actual rotor trailed tip vortices are large but finite, are concentrated at the slipstream boundary as 
“swirl,” and increase with time. Equation C3 provides an insight into observations of a large particle 
void seen in smoke flow visualization studies of rotor vortex filaments in that a finite axial flow is 
predicted for the remainder of the region r < r 0 . 
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Axial Flow Along Slipstream Boundary for t —> °° 


The angular momentum in rotor vortex filaments, as they decay, is transformed to axial 
(z coordinate) momentum. Turning to the Navier-Stokes momentum equation for v z 


3v z 3v z 1 3v 


dv. 


i ap 


+ V r — ^ + -V0^ + v z ^ = --^ + vV z v z 

at 3r r 30 3z p dz 


(C4) 


Now, as previously noted, and continuing to assume inviscid flow 


aP , 3v z 1 3v z 

— = 0 and — ^ = 

3z dz v z dt 


Therefore 


. dv, 1 dv 

2-r^- + — + V 


3v 


Z _ 


at 


ae r a r 


=o 


(C5) 


An asymptotic solution will now be derived for time approaching infinity, i.e. This 

asymptotic solution is applicable only for r > r c . 

Note that as the vortex core collapses, then 


dv z 

dt 


■0 


and 


Therefore, it can be readily seen that 


f (t) -» 0 


k = 0 as t — > 


v z — > Ae bv as t — > 


where A and b are arbitrary constants. Also, as t — > then 


\j r —> ViSgn(0 - 7t)rsin0 (as the vortex core collapses) 


Therefore, the axial velocity distribution must approach 


v z -» Ae 2 


Ae -ViSgn(0-7t)rsin0 


(C6) 
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for 


t — > oo 

The constants in the asymptotic expression, equation C6, can be approximately determined as 
follows. First, the constant b can be arbitrarily selected (for reasons that will become clear shortly) to 
be b = 2/vj. Having made that selection, then, a value for the constant A can be evaluated by noting 
that the kinetic energy of the rotor blade root vortex should be counterbalanced by the kinetic energy 
of the asymptotic filament axial velocities (equivalent to swirl velocity in rotor wake global 
coordinates versus the filament local coordinates). Therefore 

K.E. Blade Boot Swirl = K-E-Asymptotic Filament/Slipstream Swirl (C7a) 



2v 2 £2r 1 2 

t =- rdrdQ = 2 A 

(Qr) 2 +v 2 _ 



e y dydx 


(C7b) 


The expression for the blade-root swirl velocity, in the left side of equation C7b, is derived in 
reference 8 by means of the method of calculus of variation for an extended momentum theory for 
rotor induced power and thrust. The kinetic energy equivalence relationship as represented by 
equations C7a-b is only approximate as only the blade root swirl within the slipstream boundary is 
accounted for. Through the above auspicious choice of the constant b in equation C6, the kinetic 
energy integral for the filament asymptotic axial velocity has a closed form solution. Another choice 
of the constant b would not have resulted in a closed form solution. 


Ignoring higher order terms as being negligible in the integration results, the following solution for 
the constant A is derived. (Note Vj = f2R.) 
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